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Cosmology enters precision testing period with the observational experiments advancement. We 
have not arrived to decipher the origin of the dominant dark energy component of our universe, 
although we have located some of its characters. In such a situation the possibilities of a simple 
four-fermion interaction model are investigated to describe dark energy characters. We explicitly 
calculate the vacuum expectation value of the energy-momentum tensor and discuss the results 
with comparison to the equation of state for dark energy requirement. It is found that the negative 
pressure is realized at low temperature. It is an appropriate property as a dark energy candidate. 
We also investigate some cosmological tests for the dark energy via its equation of state. 
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I. INTRODUCTION 

Cosmological problems are as old as human history, 
but cosmology is a relatively new science. In recent 
years two main observational results have revealed a pro- 
found implication in our current knowledge of the uni- 
verse. These are the confirmation of the flat geometry 
and its accelerating expansion of our universe in the rel- 
atively recent epoch, i.e., with redshift of roughly 0.5 or 
less, driving by an unidentified component as coined Dark 
Energy by M. Turner [ij. 

The recent observations have dramatically changed our 
previous intuition. The evolution of the universe was 
slowing down due to gravitational attractive force and 
also alleviated the universe age crisis problem with re- 
spect to the oldest stars and globular clusters in our 
Milk Way Galaxy. Respectively, the first acoustic peak 
of the temperature fluctuations on the cosmic microwave 
background radiation (CMBR) is observed by the exper- 
iments of BOOMERANG, MAXIMA and WMAP @j. It 
favors a spatially flat universe and indicates our universe 
dominated by about 2/3 dark energy and 1/3 matter 
(with mostly the named non-relativistic cold dark mat- 
ter) components. The direct evidences from measure- 
ments of luminosity-redshift relations of SN la Q sug- 
gest that there exists a mysterious component so called 
dark energy with a negative pressure. It is as large as 
an anti-gravity contribution to accelerate the universe 
expansion. A simple candidate which shows the proper- 
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ties of the dark energy is found in the vacuum energy 
or cosmological constant. However, the energy scale to 
the cosmological constant is about 120 orders of mag- 
nitude smaller than that of the vacuum energy at the 
Planck time. Some fine tuning of parameters is neces- 
sary in almost all the models of particle physics. Thus 
a naturalness problem exists in some simple candidates. 
Furthermore, recent analysis of the observational data 
from WMAP combined with 2dF, SDSS, DASI and the 
previous CMB observations |4] disfavors the simplest cos- 
mological constant explanation. 

When describing the evolution of the universe, nor- 
mally we start with solving the Einstein equation. Since 
we know that recent observations favors the flat, homo- 
geneous and isotropic universe in Large Scale Survey, it is 
a good approximation to use the energy-momentum ten- 
sor for the dominant matter or radiation mainly depends 
on time or temperature when we consider the global 
evolution of the universe. For the perfect fluid model, 
the energy-momentum tensor consists of the total en- 
ergy density p and smoothly distributed pressure p by 
T afj = diag(p,p,p,p). 

The isotropic and homogeneous universe can be expli- 
cated by Friedmann-Robertson- Walker (FRW) metrics. 
In this case Einstein equation reduces to only two rela- 
tions. One is Friedmann equation : 
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where a(t) is the scale factor in the FRW metrics. The 
other is the acceleration equation 
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These equations are very successful to describe the evo 
lution history of the universe from the " First Three Min- 
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utes" so far as we have observed. Furthermore, these 
can accommodate the relatively recent accelerating ex- 
pansion of our universe if we accept the existence of a 
kind of vacuum energy or the cosmological constant-like 
term. On the other hand we have another solution to 
introduce an entity endowed with the sort of dark energy 
PJ. Its energy density can govern the cosmological evo- 
lution in Friedmann equations of standard cosmology [l| . 
With the present observational data we can not finally 
distinguish yet whether the accelerating expansion comes 
from the cosmological constant or not. 

The cosmological constant can be absorbed into the 
energy density p and pressure p by a re-definition in the 
above equations. To reconcile with the accelerating ex- 
pansion of the universe we must have 
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The lower bound, p/p = — 1, corresponds to the cos- 
mological constant. When we consider the dark energy 
model that causes the accelerating expansion, the dom- 
inant substance has to satisfy this relation by its own p 
and p. 

Lots of previous works are reported in concentrating 
on a dynamical scalar field such as quintessence, spin- 
essence and so on. In such investigations a slowly varying 
energy density has been employed to mimic an effective 
cosmological constant 0- As far as we know a very lit- 
tle number of works have been done in the study of the 
fermion field as the possible dark energy candidate. 

One of the interesting candidate of the dark energy 
is postulated in the model with a fermion field. In 
some strong coupling fermion models dynamical symme- 
try breaking takes place by the fermion and the anti- 
fermion condensation below a critical temperature. As 
is well known, the chiral symmetry for quarks is bro- 
ken down dynamically in QCD. The composite operator 
of the fermion and the anti-fermion develops the non- 
vanishing vacuum expectation value. It can change be- 
haviors of the energy density and pressure of the uni- 
verse at low temperature as the relevant energy scale in 
the model can be adjusted accordingly. Therefore we 
launched our plan to make a systematic study of the 
strong coupling fermion models at finite temperature for 
a dark energy candidate. 

Here we take a four-fermion interaction model as one 
of such prototype models and discuss the possibilities to 
describe the dark energy characters by using the thermal 
quantum field theory technique ■ The spacetime curva- 
ture also has an important effect on the phase structure 
of dynamical symmetry breaking in the early stage of the 
universe 0,0,0]. However, we consider the era when 
the main contribution to the evaluation of the energy- 
momentum tensor comes from the thermal effect. We 
assume that the spacetime curvature and the expansion 
of the universe do not have a large influence to the evalu- 
ation of the energy-momentum tensor. As is known, the 
four-fermion interaction is non-renormalizable. We re- 



gard the model is an effective theory which is stemming 
from a more fundamental theory at a high energy scale. 

In the present paper we study a four-fermion inter- 
action model as a dark energy candidate. By using 
the imaginary time formalism we calculate the energy- 
momentum tensor in our model at finite temperature. 
The EOS is obtained from the energy-momentum tensor. 
Behaviors of the deceleration parameter are discussed to 
see the dark energy property of the four-fermion interac- 
tion model. Throughout this paper we use natural units 
with c = ks = h = 1 and assume the universe is geomet- 
ric flat as astrophysical observation favors. 



II. ENERGY-MOMENTUM TENSOR IN A 
FOUR-FERMION INTERACTION MODEL 

First we briefly review the definition of the energy- 
momentum tensor in the quantum field theory of fermion. 
The four-fermion interaction model is one of the simplest 
model where dynamical symmetry breaking takes place 
[Til IT^I . In D-dimensional spacetime it is described by 
the action : 



S = d D x det(V) 
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where the A is a coupling constant, N is the num- 
ber of fermion flavors and V a M is a vierbein that re- 
lates the metric and the local Lorentz metric r/ a b 
by = VuVv r lab- Below we omit the flavor index k 
of fermion field. The action (JIJ has the global chiral 
symmetry under the transformation ip — > e ll&6r ij). This 
symmetry prevent the fermion from having mass term. 
Applying the auxiliary field method we rewrite the ac- 
tion as 



S = d u xdet{V) 
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The non-vanishing vacuum expectation value of a de- 
notes the existence of the fermion and anti-fermion con- 
densation, (cr) ~ —jf^ptj}) ^ 0. It is the signal for chiral 
symmetry breaking. If the chiral symmetry broken down 
dynamically, the fermion acquire a mass term propor- 
tional to (cr). 

The classical definition of the energy-momentum ten- 
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sor is |l_3| 
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Note that the labels a, b refer to the local Lorentz frame 
associated with the normal coordinates, while /i, v denote 
the general coordinate system x^. In the quantum field 
theory the expectation value of is defined in the sense 
of the path integral formalism as 



V yj Jdet[V(x)] S V^(x) 
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where Z is the generating functional which is defined by 

Z = J V([- 9 Y/H) v([-g]^) ^([-.g] 1 /^) e* s ». 

(9) 

We would like to discuss the evolution of the energy- 
momentum tensor to study a fermionic dark energy 
model. As a first step we confine ourselves to the 
Minkowski spacetime, here. The Friedmann-Robertson- 
Walker metric for the spatial flatness is conformally re- 
lated to a Minkowski line element. For the Minkowski 
metric the generating functional is give by 
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where T[a] is the effective action : 
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The energy-momentum tensor (JSJ contains a zero-point 
vacuum energy and is divergent. We use the point split- 
ting prescription to regularize the expectation value of it 
[ilf . Substituting the Eqs.lO and JTTJ into Eq.© we 
obtain 



tensor of the free fermion with mass a. Therefore the 
energy-momentum tensor of the four-fcrmion interaction 
model is given by the sum of the quadratic potential 
term of <r, i.e., the first term in the right hand side of 
Eq. (|12() . and the energy-momentum tensor of the mas- 
sive free fermion in the leading order of 1/iV-expansion. 
It should be noted that the fermion mass a is a dynamical 
variable. We can determine the dynamically generated 
fermion mass by evaluating the effective potential. 



III. ENERGY-MOMENTUM TENSOR AT 
FINITE TEMPERATURE 

Next we study the thermal effect to the energy- 
momentum tensor. As is shown in Ea. l|12[l . the vacuum 
expectation value for the energy-momentum tensor is de- 
scribed by the free fermion propagator. Following the 
imaginary time formalism, the fermion propagator at fi- 
nite temperature is obtained from the one at T = by 
the Wick rotation and the replacements 0, 0] 
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where ft = l/T. The discrete variable uj n is given by 
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according to the anti-periodic boundary condition of the 
fermion field. After these replacement and some calcula- 
tions the energy-momentum tensor 11 12 ti at finite temper- 
ature in the leading order of the 1/N expansion reads 



(T 00 ) = Na 2 /2X 
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and 
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in the leading order of 1/iV-expansion where d x denote 
the differentiation with respect to x and Sp(x, x';a) is 
the Feynman propagator for the free fermion with mass 
a which satisfies 
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The second term in the right hand side of Eq. I|12f) corre- 
sponds to the expectation value of the energy-momentum 



where trl is the trace with respect to the fermion legs in 
an arbitrary dimensions D, trl = 2 D / 2 . Since the zero- 
point of energy-momentum tensor is unknown parameter 
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FIG. 1: Behavior of the expectation value of a for A = 
as a function of T. ero is the expectation value 

at T = 0. 



in the quantum field theory, we subtract the expectation 
value for T = 0, (0|T MI/ |0), to normalize it. Since the four 
fermion interaction is not renormalizable, we consider the 
spacetime dimension less than 4 to regularize the theory. 
Here we regards the theory for D = 4 — 2e with e suf- 
ficiently small positive as a regularization of the one in 
four dimensions. There is a correspondence between the 
positive regularization parameter e and the cut-off scale 
A of the theory [Toj . 
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where p is a renormalization scale. 

To obtain the expectation value of a which corresponds 
to the dynamical fermion mass we evaluate the effective 
potential V(a) of the four- fermion interaction model at 
finite temperature. In arbitrary spacetime dimensions 
2 < D < 4 it is given by 15] 



V(a) = a 2 /2\ 
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at the large N limit. The expectation value of a is de- 
termined by observing the minimum of the effective po- 
tential. In the present paper we are interested in the in- 
fluence of the chiral symmetry breaking. Hence we take 
the coupling A larger than the critical one and numer- 
ically calculate the behavior of a as is shown in Fig. ^ 
The dynamical fermion mass is generated below a critical 
temperature, T c . Details of the generated fermion mass 
is disscussed in Ref. for 2 < D < 4. 

Inserting the expectation value of a into the Eqs. 1181) 
and (|19|l we obtain the energy-momentum tensor of the 
four-fermion interaction model at finite T. 



For the case of chiral symmetry restoration (i.e. T < 
T c ) the fermion is massless. We get the traceless energy- 
momentum tensor 

r a/9 = diag(p,p/(D-l),p/(D-l),p/(D-l),...), (22) 
where the density p is introduced by 
Ntrl 
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which is always positive by numerical calculations and 
formally satisfies the photon like equation of state at high 
temperature, w = p/p = 1/(D — 1), with D — 4 — 2e. 
The behavior for w = p/ p = 1/ (D — 1) is very reasonable 
because a = means the massless field model! So its 
equation of state must be similar to the photon's. 

For the case of chiral symmetry breaking with a 
nonzero, we have 

T a p — parjap 

+dmg(p,p/(D-l),p/(D-l),---) 
+diag( CO ( ( 7),0,0,---), (24) 

where the small correction 

. . Ntrl 
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is also positive and the energy-momentum tensor ele- 
ments also satisfy the EOS for the dark energy with the 
background of temperature dependent cosmological con- 
stant p = N[a(T = 0) 2 - cr(T) 2 ]/(2A). 

The behavior of the energy-momentum tensor is illus- 
trated in Figs. [21and[31 As is explained above, zero-point 
of (T M!y ) can not be determined in the quantum field the- 
ory. We normalize it which satisfies (T^) = at T = 0. 
As is clearly seen, the negative pressure is realized for a 
low temperature region. The lines turn sharply at the 
critical temperature T c . 

In the following section we turn to the the classical 
cosmological tests of EOS for our model. 



IV. DYNAMICAL EVOLUTION IN THE FRW 
COSMOLOGY 

In addition to the Friedmann equation we have the 
important equation to satisfy when we solve the evolution 
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FIG. 2: Behavior of the energy-momentum tensor Too for 
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of the universe. It is the conservation for the energy- 
momentum tensor of the cosmic fluid given by 

Pi + 3H{ Pi + Pi ) = 0, (26) 

where subscript i describes each species of the fluids in 
the considering universe. Here we mainly discuss the 
evolution of the universe which is dominated by the dark 
energy of our model. 

Since the recent observational data strongly suggest 
that our universe is spatially flat [lri Il7l Il8| . that is 
k = in Eq.Q), we will also investigate only this case 
here. In this framework we mainly discuss several clas- 
sical observational tests |l9j for this mysterious dark en- 
ergy problem. 

A dimensionless quantity is more useful to discuss 
physical phenomena of our model. In the standard cos- 
mology, for convenience, the density fraction of the world 
is expressed in terms of a critical density which is defined 
by the total energy density, i.e. p c — 3Hq/8ttG. When 
we treat the flat, homogeneous and isotropic universe, 
the comoving coordinates is convenient to describe the 
dynamical evolution of it. In this treatment the coordi- 
nates at any place do not vary with the evolution of the 
universe. It is described by a scale factor a(t) that is 
a function of only cosmic time in the homogeneous and 



isotropic universe. If we find the relationship between the 
scale factor and the cosmic time, then the evolution of the 
universe is clearly understood. For example, we can di- 
rectly deduce the Hubble parameter and the deceleration 
parameters as a function of time, the dependence of the 
luminosity distance on redshift and so on. By investigat- 
ing its dynamical evolution in the dark energy dominant 
era, we compare the observable quantities derived in our 
model with the experimental data. 

We first consider a general case where equation of state 
is written by 

P = wp, (27) 

The parameter w depend on the spacetime dimension D 
and this dependence is assumed to be extremely small fol- 
low temperature region. By using the continuity equation 
(|26|l and the Ea. (|27|) . we easily obtain the evolution of 
the energy density p which obeys 

p = Po a- n , (28) 

where po is the present energy density of the universe and 
n is given by n = 3(1 + w). Substituting Eq.J2BJ in the 
Friedmann equation, we get 

H 2 = Hi [n t a- n + fl M a- D+1 + n R a- D ] , (29) 

where Hq is the Hubble parameter today and f2 a is the 
density parameter defined by fl a = p a o/pc- Since we 
consider the flat universe, the density parameters must 
satisfy the relation, J^i &i = 1- 

It is hard to solve the Friedmann equation l|29() in gen- 
eral situations, because we do not know the EOS for the 
dark energy. From Ea. H29l) . however, we can easily ob- 
tain the evolution of the scale factor. If we consider the 
dark energy dominant situation, the scale factor behaves 
as 

a(t) = ^MK t y\ for n>0 

a(t) cx exp (^jH%n t ?j . for n = (30) 

The latter solution shows the same behavior with the in- 
flationally expansion when the energy scale is sufficiently 
high. We may use it to describe the accelerating expan- 
sion with a cosmological constant in the low energy scale, 
as observed today. 

The dynamical evolution factor is just the function of 
EOS which is shown in Fig. 0] Comparing the lines in 
the figure, we can see that for the smaller value of n the 
scale factor a(t) increases faster. 

V. THE DECELERATION PARAMETER 

Following the fact that our universe is expanding, one 
may ask a question that the expansion of our universe is 
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accelerating or decelerating now, or how fast it expands? 
A parameter to describe the evolution of the universe is 
not only the Hubble parameter but also another impor- 
tant parameter, so called the deceleration parameter q 
which is defined by 



The parameter q is dimensionless parameter which de- 
scribes the rate for slowing the expansion. By using 
Eq.(j2nj) we obtain 

g=i(n-2) = i(3«; + l). (32) 

As a first consideration we assume that the state equation 
parameter w is constant in the dark energy dominant era. 

In the four-fermion interaction model we show the evo- 
lution of the deceleration parameter in Fig. |3J As is 
clearly seen in Fig^w is not a constant. The evolutions 
of w depend on both the temperature and the regular- 
ization parameter e which is related to the cut-off scale 
of the theory. However, we think that the constant w 
is a good approximation when we consider the present 
evolution of our universe because our universe is bathed 
in enough low temperature now. Generally speaking, the 
deceleration parameter q is a function of both the present 
density flj and w. Furthermore, it is also a function of the 
parameter e because we consider a non-renormalizable 
model as a low energy effective theory. Figure |3] shows 
that the relationship between the parameter w and the 
regularization parameter e. The present w (low temper- 
ature region) seems almost a monotone function of the 
parameter. The parameter w smoothly falls down and 
the curves become flat as T decreases. At the lower tem- 
perature limit w approaches the cosmological constant 
limit. It shows the exponential expansion of the present 
universe. As is shown in Eq. I|22|l w becomes l/(D — 1) 
at the high temperature limit. It is a property of the 
massless fermion. Therefore the exponential expansion 
only appears at the low temperature region in our model 
without any fine tunings. It is a very important property 



to a dark energy candidate. Slopes of each lines in Fig[S] 
are depend on the coupling A. For larger coupling A the 
slope becomes steeper. However, the behaviors of w does 
not change qualitatively as A varies. 

From the recent SN la and Cosmic Microwave Back- 
ground observations [j| a cosmological constant-like dark 
energy may exist today which accelerate the expansion of 
the universe. It suggests that the decelerating parameter 
lies in a range, — 1 < q < 0. This observational con- 
straint is naturally satisfied in the four-fermion interac- 
tion model. Recently there are many dark energy models 
of simple scalar fields like the quintessence |6(, however, 
they all need fine tuning of some parameters to certain 
extents. In the present model, on the other hand, fine 
tuning is not necessary. In other words the model has 
some properties which are realized in the quintessence 
scenario without any cosmological constant-like terms. 
Certainly we need more SN la observational data to con- 
firm the important constraint, which is just the main 
mission of the near future SNAP project. 



VI. CONCLUDING REMARKS 

We investigated the cosmological effects of the vacuum 
expectation value of energy momentum tensor in a four 
fermion interaction model. According to the imaginary 
time formalism we derive the explicit form of the en- 
ergy momentum tensor. The energy-momentum tensor 
depends on the temperature and the dynamically gener- 
ated fermion mass. The fermion field is massless and has 
the similar property as the radiation field above the criti- 
cal temperature. We consider that the temperature of the 
present universe is much less than the critical one and the 
chiral symmetry is broken down dynamically. Thus the 
composite operator of the fermion and the anti-fermion 
develops a non-vanishing vacuum expectation value in 
the present stage of the universe. The vacuum expecta- 
tion value plays a role of the cosmological constant-like. 

Applying it to models of the universe without an ex- 
plicit cosmological constant term, we discussed the de- 
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pendence of the scale factor a(t) on the state equation 
parameter w and evaluate the deceleration parameter q. 
We found that the behavior of the EOS is dramatically 
changed at low temperature. A dark energy property 
was naturally obtained in the four-fermion interaction 
model. In the present analysis the fcrmion mass scale can 
be at the scale of neutrino, QCD, electroweak symmetry 
breaking, extra-dimension [20| and so on. We have had a 
loosely constrained parameter range for the cosmological 
constant [2]]]. For more concrete analysis we must apply 
the results to some cosmological tests with comparisons 
to observational data, which are under our investigations. 

Our model is only a prototype one, but we believe that 
the results are worth discussing further about the real 
world. By observational data analysis we know the dark 
energy is almost smoothly distributed 22] , which implies 
a kind of very weak interaction that we are not very clear 
about yet. Here we just employ the simple four-fermion 
interaction model to give the global energy momentum 
tensor which depends on temperature and some model 
parameters. It might be too simple, however we found 
interesting results. As a next step, it is very important 
to apply our analysis to other well-defined and more real- 
istic models. One of the interesting models is SUSY NJL 
model HHiil. For the SUSY NJL model the effective 
action is written as 

r[cr] = J d 4 x[Tr\n(i - a) 

-2Trln(<9 2 - a 2 ) + Na 2 /(2X)] . (33) 

According to the supersymmetry, contributions from the 
fermion and boson loops cancel out each other exactly 
at T = 0. Since the thermal effect restore the chiral 



symmetry, dynamical symmetry breaking does not take 
place even at finite temperature. If we introduce a soft 
SUSY breaking term such as A(f>*(j> and/or A(p c *(j> c in the 
action 25], then the contribution from the boson loop 
is modified. The composite operator ipip can develop a 
non-vanishing vacuum expectation value which satisfies 
the dark energy property to some range. The scalar mass 
depends on the coupling parameter. When the coupling 
parameter is large enough, A > 2er, the scalar mass m = 
s/ (cr 2 — A) becomes heavier than the dynamical fermion 
mass. We can choose the coupling parameter A to realize 
the appropriate scale for the cosmological constant and 
the negative pressure of the present universe accordingly. 

So far the origin of the dark energy in the current 
cosmological models is still under discussions with many 
theoretical speculations. It is accommodated in the stan- 
dard Friedmann cosmology based on the relativistic grav- 
itational theory with a cosmological constant-like term. 
It can also be interpreted in extended gravity scenarios 
and corresponding cosmological models, like the brane 
cosmologys. That attracts us to further research works 
continuously. 
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